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ABSTRACT 

The Yale stellar evolution code has been modified to include the combined effects 
of diffusion and rotational mixing on ^H, "^He and the trace elements "^He, ^Li, ''Li, and 
^Be. The interaction between rotational mixing and diffusion is studied by calculating 
a number of calibrated solar models. The rotational mixing inhibits the diffusion in the 
outer parts of the models, leading to a decrease in the envelope diffusion by 25 - 50%. 
Conversely, diffusion leads to gradients in mean molecular weight which can inhibit the 
rotational mixing. The degree to which gradients in mean molecular weight inhibits the 
rotational mixing is somewhat uncertain. A comparison to the observed solar oblateness 
suggests that gradients in the mean molecular weight play a smaller role in inhibiting the 
rotational mixing previously believed. This is reinforced by the fact that the model with 
the standard value for the inhibiting effect of mean molecular weight on the rotational 
mixing depletes no Li on the main sequence. This is in clear contrast to the observations. 
A reduction in the inhibiting effect of mean molecular weight gradients by a factor of 
ten leads to noticeable main sequence Li depletion. 

Subject headings: stars: interiors - stars: rotation - stars: abundances - sun: rotation 
- sun: interior 
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1. Introduction 

In standard stellar models, it is assumed that no 
material is transported in stellar radiative regions. 
^Li is a particularly good tracer of particle trans- 
port in stellar radiative zones, because it is destroyed 
at a temperature of ~ 2.5 x 10® K. For stars with 
M ^ 1.0 Mq, the Li destruction zone is typically lo- 
cated in the radiative region of the star throughout 
the pre-main sequence and main sequence evolution. 
Thus, if no particle transport occurs, then the sur- 
face value of Li that we observe should be the initial 
amount. The solar photospheric Li value is a factor 
of 200 smaller than that observed in meteorites (An- 
ders & Grevesse 1989), indicating that some particle 
transport has occurred in the solar radiative regions. 
A comparison of Li abundances in the Hyades with 
younger clusters also provides strong observational 
evidence for main sequence Li depletion (Thorburn 
ei al. 1993; Soderblom ei al. 1993). 

The most striking evidence of particle transport in 
stellar radiative regions is found in stars in the tem- 
perature range 6400 - 6800 K. These stars are defi- 
cient in Li by 1 - 2 orders of magnitude (Boesgaard & 
Tripicco 1986, Boesgaard & Budge 1988). Other clus- 
ters of similar age, and older, also contain a Li gap: 
Coma, UMa, and NGC 752. However, the younger 
clusters Pleiades and a Per do not show the pres- 
ence of a Li dip (see discussion by Boesgaard 1991). 
To date, no stellar evolution models have been able 
to fully explain the presence of the Li gap, though 
attempts have been made (Pinsonneault, Kawaler & 
Demarque 1989; Richer & Michaud 1993). Its exis- 
tence clearly points to rapid transport in the radiative 
regions of some stars. 

Another difficulty for the standard stellar evolution 
models was found by recent observations of Pleiades 
and Hyades stars which demonstrate that a large dis- 
persion exists in the Li abundance for stars of the 
same mass (Soderblom et al. 1993; Thorburn et al. 
1993). In standard stellar models, the amount of 
^Li depletion depends only on the star's mass, compo- 
sition and age. Thus, another parameter must govern 
the depletion of Li. 

Together, these observations require that some 
transport mechanism is operating in stellar radiative 
regions. In this paper, we examine the hypothesis 
that this particle transport is caused by microscopic 
diffusion and rotational instabilities and explore the 
interactions between the two processes. Previous 
studies have examined particle transport due to dif- 
fusion (eg. Michaud 1986; Proffitt k Michaud 1991; 
Chaboyer et al. 1992; Richer & Michaud 1993), due to 



rotational instabilities (Endal & Sofia 1978; Pinson- 
neault et al. 1989, 1990), and due to the combined 
effect of diffusion and advection due to meridonal cir- 
culation (Charbonneau & Michaud 1991; Charbon- 
neau 1992). In this paper, we study the combined ef- 
fect of diffusion, and rotation induced transport due 
to: meridonal circulation, dynamical shear instabil- 
ities, the Solberg-Hoiland instability, the Goldreich- 
Schubert-Fricke instability and the secular shear in- 
stability. In contrast to Charbonneau & Michaud, 
who treated the meridonal circulation as a advec- 
tion process, we treat meridonal circulation as a dif- 
fusion process. This is in agreement with the work 
of Chaboyer & Zahn (1992), who showed that appro- 
priate rescaling of the turbulent diffusion coefficient 
allows the combined advection/diffusion process in- 
duced by meridonal circulation to be treated solely 
as a diffusion process. 

It is important to note that diffusion will not just 
affect the surface abundance of ^Li, for diffusion will 
also cause '^He to sink relative to hydrogen. He- 
lium diffusion is particularly important in affecting 
the structure (Vauclair, Vauclair & Pamjatnikh 1974; 
Bahcall & Pinsonneault 1992; Richer, Michaud & 
Proffitt 1992) and evolution (Noerdlinger & Arigo 
1980; Chaboyer et al. 1992) of a star because he- 
lium is very abundant and its diffusion time scale is 
relatively short. 

The rotational models at Yale have been able to 
reproduce many of the observations of ^Li depletion 
patterns and rotational velocities observed in young 
cluster stars (Pinsonneault et al. 1990). However, 
these models ignored microscopic diffusion, which also 
can have important effects on the calculated particle 
transport. In particular, the diffusion of '^He will lead 
to gradients in the mean molecular weight. A gra- 
dient in mean molecular weight tends to inhibit ro- 
tational mixing, as extra energy is required to trans- 
port material to regions of different mean molecular 
weight (Mestel 1953). Thus, diffusion will tend to in- 
hibit the mixing induced by rotation. Conversely, the 
rotational mixing can inhibit the separation process 
caused by microscopic diffusion or radiative levita- 
tion. When discussing the competing effects of dif- 
fusion and rotational mixing, one must keep in mind 
the key difference in these effects: diffusion is a sepa- 
ration process, with elements sinking with respect to 
hydrogen; while rotational mixing does not separate 
elements, it merely mixes them into different regions. 

The Sun provides a unique calibration and test 
of stellar models since its properties (age, chemical 
composition, luminosity, etc) are known far more ac- 
curately than any other star. In addition, helioseis- 
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mology and the solar oblateness provide information 
regarding the internal rotation of the Sun which is 
not available for other stars. Finally, by comparing 
photospheric and meteoritic abundances it is possible 
to determine the amount of destruction of ^Li. The 
Sun is the only star for which it is possible to make 
this measurement. For this reason, we start our ex- 
ploration of the effects of microscopic diffusion and 
rotational mixing on stellar evolution models by con- 
structing a series of calibrated solar models. It must 
be stressed however, that in performing this research, 
we have also considered the constraints from observa- 
tions of Li abundances and surface rotation velocities 
in open cluster stars. Observations of open cluster 
stars of various ages and masses provide provide cru- 
cial information regarding the time-scale and mass 
dependance of ^Li depletion which are not available 
solely from a study of the Sun. 

For comparison purposes a number of solar models 
have been constructed: standard solar models; mod- 
els which only includes diffusion; a model which only 
includes rotational mixing; and models which include 
both rotational mixing and diffusion. This paper is 
the first in a series which will examine the role of par- 
ticle transport in stellar radiative regions. Here, we 
introduce the method which is used to model the par- 
ticle transport, and concentrate on the interaction be- 
tween diffusion and rotational mixing by studying the 
calibrated solar models. Paper II (Chaboyer, Demar- 
que & Pinsonneault 1994) explores the various param- 
eters associated with the stellar models and compares 
the models to observations of Li and rotation veloci- 
ties in young cluster stars. A related paper (Chaboyer 
& Demarque 1994) examines the constraints that ob- 
served halo star Li abundances put on particle trans- 
port in stellar radiative regions. 

In §2. we provide a description of the construc- 
tion of our solar models, including the treatment of 
particle in the radiative regions of our models, and 
list the observed solar parameters which are used to 
calibrate our models. The calibrated solar models 
are presented in §3. and the interaction of rotational 
mixing and microscopic diffusion is discussed. The 
^Li depletion and surface rotation history of the solar 
models are shown in §4. Finally, the main results of 
this paper are summarized in §5. 

2. Model Construction 
2.1. Overall Approach 

The solar models are constructed by evolving a 
fully convective, pre-main sequence 1 Mq model (with 
a nominal age of 1 Myr) to the age of the Sun. The 



initial X , Z and a are varied until the solar age model 
has the solar values of luminosity, radius and Z/X. 
The one dimensional rotating stellar structure equa- 
tions (see Pinsonneault ei al. 1989) are solved using 
the Henyey relaxation method and are converged to 
better than 0.1% in each time step. Typically, ~ 330 
evolutionary time steps (of which 70 were on the main 
sequence) were required to evolve to a solar model. 

After each evolutionary time step, the stability of 
the angular velocity profile from the previous time 
step is checked and the diffusion coefficients are cal- 
culated (see §2.3.). If the diffusion coefficients are 
large, then a series of small time steps are taken, with 
the thermal structure determined by linearly interpo- 
lating between the previous and current model. In 
each time step, the coupled nonlinear diffusion equa- 
tions for angular momentum and chemical composi- 
tion (equations (1) and (3)) are solved using an iter- 
ative technique, with the convergence criterion set to 
0.1%. 

All of the solar models presented here match the 
solar radius and luminosity to within 0.1%, while the 
surface Z/X was matched to the observed value to 
within 1.0%. Models which include rotation were re- 
quired to match the observed solar rotation rate and 
^Li depletion to within 1.0%. This was achieved by 
varying the values of and (see equations (3) - 
(5)). 

Our models include the latest physics available to 
us. The energy producing nuclear reaction rates are 
from Bahcall & Pinsonneault (1992); reaction rates 
for the ®Li, ^Li and ^Be are from Caughlan & Fowler 
(1988); the high temperature opacities from Iglesias & 
Rogers (1991); the low temperature opacities (below 
lO'* K) are from Kurucz (1991, and private commu- 
nication); while the surface boundary conditions are 
determined using Kurucz model atmospheres (Kurucz 
1992, and private communication; Howard, 1993). 
For temperatures above 10® K, a relativistic degen- 
erate, fully ionized equation of state is use. Below 
10® K, the single ionization of ^H , the first ionization 
of the metals and both ionizations of '^He are taken 
into account via the Saha equation. 

2.2. Calibrated Solar Models 

A calibrated solar model is a 1 Mq stellar model 
which, at the age of the Sun, has the observed solar 
radius, luminosity and surface chemical composition. 
The age of the Sun may be determined by radioac- 
tive dating of the oldest meteorites (Tilton 1988), 
modern theories of star formation, observations of T- 
Tauri star and pre-main sequence calculations. The 
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best estimate for the age of the Sun is 4.52 ± 0.4 Gyr 
(Guenther et al. 1992, Guenther 1989). We have 
adopted 4.55 Gyr for the solar age in all of our so- 
lar models. For the radius, mass and luminosity, we 
adopt the same values as Guenther et al. (1992): 
Mq = (1.9891 ± 0.0004) X 10^3 g; Rq = 6.9598 x 10^° 
cm at r = 2/3; and Lq = 3.8515 x lO^^ergs s'^. 

The chemical composition of the Sun may be in- 
ferred by examining the photospheric abundances in 
the present Sun and the abundances found in the old- 
est meteoritics (Anders & Grevesse 1989). For most 
elements, there is excellent agreement between the 
two techniques. The exceptions are ^Li, ^Be and Fe. 
It now appears that the Fe abundances in the photo- 
sphere were incorrect; new measurements give similar 
abundances to those found in meteorites (Biemont ei 
al. 1991; Holweger ei al. 1991). In our solar models, 
we adopt the Anders & Grevesse (1989) photospheric 
heavy element mixture with the meteoritic Fe abun- 
dance. 

Unfortunately, it is impossible to determine the 
abundance of or '^He, and so the mass fraction 
of (X) is treated as a free parameter. The mass 
fraction of '^He (Y) is determined by requiring the to- 
tal mass fraction (X + Y + Z , where Z is the mass 
fraction of the heavy element mixture) be equal to 
one. The Anders & Grevesse photospheric mixture 
with meteoritic Fe has Z/X = 0.0267 ± 0.001. 

There is one other free parameter in a standard 
solar model, the mixing length, a which is used to 
describe the convective transport of energy. The out- 
ermost layers of the Sun are convective, thus the effi- 
ciency of the convective energy transport largely de- 
termines the radius of our solar models. The value of 
a is determined by the requirement that the models 
have the observed solar radius at the solar age. 

In addition to the above surface properties of the 
Sun, it has recently become possible to determine 
the location of the base of the convection zone us- 
ing p-mode oscillation frequencies. The depth of 
the convection zone is found to be 0.713 ± O.OO3i?0 
(Christensen-Dalsgaard, Gough & Thompson 1991). 
For a given set of input physics (opacities, atmo- 
spheres, diffusion, etc), the depth of the convection 
zone is a fixed quantity and can be considered a pre- 
diction of that particular solar model. 

For the rotating models, the initial angular mo- 
mentum distribution must be prescribed. Since the 
models are evolved starting with a fully convective 
pre-main sequence model and we assume solid body 
rotation within convection zones^ , the initial angu- 

^ The validity of this assumption is addressed at the end of this 



lar momentum distribution is fully described by the 
initial surface rotation velocity. Studies of T-Tauri in- 
dicates that 1 Mq stars have rotation periods in the 
range 1.2-12 days and rotation velocities between 8 - 
70 km/s (Bouvier 1991, Bouvier ei al. 1993). The ma- 
jority of T-Tauri stars have rotation velocities around 
10 km/s. It appears that the Sun is over-depleted 
in ^Li compared to other stars of its age. Since the 
amount of mixing in a star increases as its initial ro- 
tational velocity increases, it is likely that the Sun 
was initially a rapid rotator. For this reason, we have 
chosen an initial rotation velocity of 30 or 50 km/s 
for the solar models. 

Two additional constraints are imposed on the ro- 
tating models. We require that the models have the 
solar equatorial velocity rate 2.02 ± 0.04 km/s (Lib- 
brecht & Morrow 1991) and the observed depletion of 
^Li. For convenience, we have calibrated the surface 
rotation velocity of our models to within 1% of 2.0 
km/s. 

The depletion of ^Li is inferred by comparing the 
cosmic abundance to the photospheric abundance and 
is a factor of 200 (eg. Duncan 1981). A comparison 
of the meteoric abundance to the photospheric abun- 
dance shows the depletion to be I40I30 (Anders & 
Grevesse 1989). For most of this work, we have as- 
sumed a solar ^Li depletion of 200. It is important to 
realize that in standard solar model, the ^Li deple- 
tion is very small (factor of 2 - 5). Thus, we do not 
attempt to fit the ^Li depletion with our non-rotating 
models. 

The Sun is unique among stars in that it is possi- 
ble to obtain information regarding the interior ro- 
tation of the Sun. The oblateness of the Sun (e) 
is a function of the surface solar rotation rate and 
of the quadrupole moment of the Sun. Since the 
quadrupole moment is a function of the interior ro- 
tation velocity, the oblateness of the Sun constrains 
the allowed rotation profiles of a solar model. Ob- 
servations of the oblateness of the Sun are rather 
difficult. Dicke & Goldenberg (1967) found e = 
(5.0 ± 0.7) X 10~®, a value which implies that general 
relativity incorrectly predicts the perihelion advance 
of Mercury. However, Hill & Stebbins (1975) deter- 
mined e = (9.6 ± 6.5) X 10~®, and suggested that the 
Dicke & Goldenberg observations were hindered by 
brightness variations between the pole and equator 
of the Sun. Dicke, Kuhn & Libbrecht (1985) found 
e = (2.0 ± 0.4) X 10"^ Most recently, Sofia, Heaps & 
Twigg (1994) measured the oblateness using a balloon 
borne experiment and found e = (8.6 ± 0.9) x 10~®. 
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This is by far the most accurate measurement of the 
solar oblateness to date. We have chosen to ignore 
the Dicke & Goldenberg (1967) results as they are in- 
compatible with the more recent determinations. In 
using the solar oblateness measurements to test the 
models, it must be remembered that the oblateness 
yields information about the gross properties of the 
rotation profile; different rotation laws may have sim- 
ilar values for the solar oblateness. 

In principle, the splitting of the solar p-modes con- 
tain detailed information on the latitude and depth 
dependance of the solar angular velocity. Libbrecht 
(1989) obtained observations of the I = 10 - 60 p- 
modes, which yield information about the rotation 
velocity down to ~ 0.4 Rq (Demarque & Guenther 
1988). These observations show that in the equato- 
rial plane, the convection zone essentially rotates as 
a solid body. As is well known, there is a signifi- 
cant dependence of the surface solar rotation rate on 
latitude. The latitude dependence of the angular ve- 
locity decreases below the surface such that at the 
bottom the convection zone the angular velocity is 
independent of latitude (Libbrecht & Morrow 1991) 
to within ~ 10%. Measurements of low I p-modes 
from space (Toutain & Frohlich 1992) yield informa- 
tion regarding the rotation rate of the deep interior 
of the Sun and suggest that at r ~ 0.2 Rq, the Sun 
rotates ~ 5 times faster than the surface (Toutain 
k FrShlich 1992; Goode, FrShlich k Toutain 1992). 
The IRIS network observed the 1=1 mode (Loudagh 
et al. 1993) and found a smaller rotational splitting 
leading to central rotation rates of ~ 3 times the sur- 
face value. We must caution however, that the weight 
of the core (r < 0.3 Rq) on the total splitting is less 
than 10%, even these low I models (Loudagh et al. 
1993). In addition, these results are determined via 
an inversion procedure, which is sensitive to obser- 
vational errors and the type of inversion procedure 
employed (see Goode et al. 1991). 

2.3. Treatment of Mixing in Radiative Re- 
gions 

In the radiative regions of a star the transport 
of angular momentum and chemical composition are 
treated in a similar manner to Pinsonneault et al. 
(1989) using the coupled nonlinear equations: 
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where uj is the angular velocity, Xi is the mass frac- 
tion of chemical species i, I/M is the moment of in- 
ertia per unit mass, Dm,i and Dm,2 are derived from 
the microscopic diffusion coefficients and multiplied 
by the adjustable parameter /„, Drot is the diffusion 
coefficient due to rotational instabilities and fc is an 
adjustable parameter, which represents the fact that 
the transport of angular momentum and a chemical 
species occur on different time-scales. The value of 
fc is approximately 0.033 (Chaboyer k Zahn 1992). 
The value of fc is varied (within reasonable limits) un- 
til the correct ^Li depletion for the Sun is obtained. 
The transport of angular momentum due to the mi- 
croscopic diffusion is not included, as it is a second 
order effect. 

In all of the models, a small over-mixing zone (0.02 
to 0.05 pressure scale heights) below the base of the 
surface convection zone has been included. This over- 
mixing zone has been included as many of the esti- 
mates for the rotational mixing coefficients go to in- 
finity at the base of the convection zone. The size 
of the over-mixing zone at the base of the convection 
zone can be important in determining the amount of 
pre-main sequence ^Li depletion. However, the rela- 
tively small over-mixing zones used in these models 
do not appreciable effect the models. 

2.3.1. Microscopic Diffusion Coefficients 

The microscopic diffusion coefficients of Michaud 
k Proffitt (1993) are used throughout this paper. 
For and '^He, these coefficients are ~ 20% larger 
than the Bahcall k Loeb (1990) coefficients. Recent 
work by Thoul, Bahcall k Loeb (1994) who solve the 
full system of Burgers equations find agreement with 
Michaud k Proffitt to within <^ 15%. We have intro- 
duced an adjustable parameter which multiplies 
all of the microscopic diffusion coefficients. For the 
work presented here, we have set = 0.8 or 1.0. 

This work does not include radiative levitation. 
This is a good approximation when the atoms are 
fully ionized. The light elements ^He, ^Li, and 
^Be are fully ionized throughout most of the ra- 
diative region of our models and so we use the above 
equations for the microscopic diffusion coefficients for 
the light elements. It is not clear how important ra- 
diative levitation (or diffusion) is for heavier elements. 
Observations of heavy elements in Hyades stars show 
similar abundances in stars of different spectral types 
(Boesgaard k Budge 1988). Thus, we have chosen 
not to include the effects of microscopic diffusion on 
any element heavier than ^Be . 
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2.3.2. Rotational Diffusion Coefficients 

Rotational mixing is triggered by various hydro- 
dynamic instabilities that occur in a rotating star 
and lead to the bulk transport of matter in the ra- 
diative regions of a star (see Zahn 1993 for a recent 
review). The evolution code at Yale is able to follow 
the evolution of rotating stars and includes the ef- 
fects of angular momentum loss due to stellar winds, 
and the transport of angular momentum and chemi- 
cal species due to various hydrodynamic instabilities 
(Pinsonneault et al. 1989). It is assumed that con- 
vection enforces solid body rotation in the convective 
regions of a star, so that the rotational instabilities 
are only effective in the radiative regions. 

Rotation induced mixing occurs on two different 
time-scales: dynamical and secular. The dynami- 
cal instabilities have associated time-scales for mixing 
which are orders of magnitude smaller than other time 
scales involved in the hydrostatic evolution of a star. 
As such, if a region is dynamically unstable, the an- 
gular momentum distribution is instantaneously read- 
justed to a point of marginal stability. If a composi- 
tion gradient existed in the unstable region, then it is 
reduced by the same amount as the angular velocity 
gradient. Angular momentum and particle numbers 
are conserved within the regions subject to the insta- 
bilities. In its present form, our stellar evolution code 
treats two types of rotation-related dynamical insta- 
bilities: the dynamical shear instabilities Zahn (1974) 
and the Solberg-Hoiland instability (cf. Wasiutynski 
1946). The stability criterion used are the same as 
Pinsonneault et al. (1989). 

Rotation induced instabilities whose time-scale for 
mixing is comparable to evolutionary time-scales (sec- 
ular), are treated by solving the coupled diffusion 
equations (1) and (3). The turbulent diffusion coeffi- 
cient Drot in these equations is taken to be the prod- 
uct of a velocity estimate v and path length r (Zahn 
1993). Note that Pinsonneault et al. (1989, 1990) 
determined the turbulent diffusion coefficient as the 
product of a velocity estimate and a local scale height. 
Thus, the estimate for the rotational mixing in these 
models is somewhat different from Pinsonneault et al. 
(1989, 1990). The secular instabilities treated in this 
work are: meridional circulation (von Zeipel 1924; 
Eddington 1925; Kippenhahn & Mollenhof 1974), the 
Goldreich-Schubert-Fricke (GSF) instability (Goldre- 
ich & Schubert 1967; Fricke 1968) and the secular 
shear instability (Zahn 1974). In all cases, the cri- 
terion for stability and estimate of the diffusion ve- 
locity are calculated in the same manner as Pinson- 
neault et al. (1991). Gradients in the mean molecular 



weight {ji) inhibit mixing due to meridional circula- 
tion (Mestel 1953) and the GSF instability. This ef- 
fect is taken into account by calculating a fictitious jj, 
velocity (Mestel 1953) which is subtracted from the 
velocity estimate for the meridional circulation and 
GSF circulation. 

The amount of mixing induced by the GSF in- 
stability is particularly uncertain. For this reason, 
we have introduced an adjustable parameter, fasF 
which multiplies the velocity estimate of James & 
Kahn (1970, 1971). Values oi fasF greater than one, 
imply very efficient mixing due to the GSF instability. 

In addition, the inhibition of rotational mixing due 
to gradients in the mean molecular weight is not well 
understood. As the microscopic diffusion builds up 
gradients in the mean molecular weight, this is a key 
parameter when studying the interaction between dif- 
fusion and rotational mixing. For this reason, we in- 
clude a parameter which multiplies the velocity es- 
timate for the fictitious jj, current and determines the 
efficiency of rotational mixing in areas which have a 
mean molecular weight gradient. Hence, the estimate 
for the particle transport caused by rotational insta- 
bilities contains three parameters: /c, fasF and 
The parameter is determined by requiring that the 
^Li depletion in our models match the observed de- 
pletion. The parameters and fasF are treated as 
free parameters in this work. 

We have not included angular momentum trans- 
port by magnetic fields or internal gravity waves 
(Press 1981). Either could in principle be important, 
but both suffer from significant uncertainties in quan- 
tifying their impact. Torsional oscillations act to re- 
move angular velocity gradients in a relatively short 
timescale (see for example Spruit 1987) provided that 
the magnetic field strengths are allowed to grow ar- 
bitrarily large. However, magnetic instabilities could 
prevent the growth of large internal fields (see Mes- 
tel and Weiss 1987), or internal magnetic fields could 
be decoupled from the envelope (see Charbonneau & 
MacGregor 1993 for a discussion). We note that both 
stellar and solar observations provide evidence for dif- 
ferential rotation with depth inside stars, which con- 
strains the role of magnetic fields and waves. Further 
data on internal stellar rotation, and further advances 
in theory, will be needed to quantify the role of these 
potentially important players. 

2.3.3. Angular Momentum Loss 

It is clear from observations of young clusters, and 
from the fact that the Sun is a slow rotator, that stars 
lose angular momentum as they evolve (Kraft 1970). 
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It is likely that this loss is caused by a magnetic stellar 
wind (Schatzman 1962). Kawaler (1988), following 
Mestel (1984) developed a general parameterization 
for the loss of angular momentum due to magnetic 
stellar winds. The rate of angular momentum loss 
was proportional to uj^+^^l^ ^ where # is a measure 
of the magnetic field geometry and may vary from 
to 2. Since Kawaler's work, it has become clear 
that such a simple dependance on to is not valid; in 
particular it appears that there is a saturation level 
in the angular momentum loss (Stauffer & Hartmann 
1987; Bouvier 1991) in the pre-main sequence. For 
this reason we have modified Kawaler's relation to be 

/ . X 1-2N/3 

/ . \ 1-2N/3 
/ M \ 4N/3 

where Wcrit introduces a saturation level into the an- 
gular momentum loss law, A'„, = 2.036x 10^^-(1.452x 
10®)^ in cgs units, M is the mass loss rate in units of 
lO"-'^'* Mq yr""'^. The dependance on mass loss rate is 
very weak, for this study we have set it to 2.0 x 10""'^ ■* 
Mq yr""'^. Thus, there are 3 adjustable parameters in 
the wind loss law, Wcrit, N and /j,. The parameter N 
is between and 2. The value of N determines the 
rate of angular momentum loss with time. The value 
of fk determines the total amount of the angular mo- 
mentum loss and it is adjusted until our models have 
the solar rotation velocity at the solar age. There are 
no reliable values for Wcrit! it is a free parameter in 
our models. 

To be consistent, a global outward drift term 
should be added to equation (3) in order to take into 
account the effect of mass loss on element diffusion. 
However, we have run a few test cases which indi- 
cate that, at the level of mass loss considered here 
(2.0 X lO""'^'* Mq yr""'^), the mass loss term has little 
effect on the element diffusion for the masses below 
1.3 Mq). For this reason, the mass loss term in the 
element diffusion equation (3) has been ignored. 



3. Solar Models 

3.1. Basic Parameters 

In order to study the interaction of rotation and 
diffusion, it is important to separate out the various 
effects. For this reason, we have constructed a number 
of calibrated solar models: standard solar models (no 
rotation or diffusion), a purely rotational model; mod- 
els which only include diffusion; and a number mod- 
els which include rotation and diffusion. Although 
there are a number of adjustable parameters associ- 
ated with the rotating models. Extensive testing in- 
dicates that only plays a key role in determining 
the interaction between rotational mixing (Chaboyer 
1993) and diffusion. 

The properties of our solar models are shown in 
Table 1. Each solar model is designated by a two let- 
ter code, given in the first column, which will be used 
to refer to the models throughout the rest of this pa- 
per. For each model, the input mixing length, helium 
abundance (Yo), Zg/ Xg and extent of the overshoot 
region are shown in columns 2-5. Column 6 indi- 
cates the constant which multiplies the diffusion coef- 
ficients. A value of indicates that diffusion was not 
included. Columns 7-11 give the rotation parame- 
ters. A value of in column 7 (the initial rotation ve- 
locity) indicates that rotation is not included in that 
model. Some model results are given in the last two 
columns. The ratio of the initial ^Li abundance to the 
final surface value is shown in column 12. Finally, col- 
umn 13 gives the convection zone depth. The different 
types of models (rotation and diffusion, pure rotation, 
pure diffusion or standard) are separated from each 
other by a horizontal line. Discussion of some of the 
models in Table 1 is deferred to Paper II. 

It is clear from Table 1 that the depth of the con- 
vection zone matches the observed value (0.713 ± 
0.003) if we include rotation and microscopic diffu- 
sion. Pure diffusion models have slightly deeper con- 
vection zones, and the match to observations is not as 
good. Models which do not include diffusion (with or 
without rotation) fail to match the observed convec- 
tion zone depth, and so might be ruled out as good 
solar models. We caution, however, that solar mod- 
els that have no diffusion but have opacities enhanced 
near the base of the convection zone can match the 
observed depth of the convection zone (Guzik & Cox 
1993). 

The present solar rotation curve is a strong func- 
tion of Figure 1 plots the solar rotation curve for 
the four rotating models which differ in All of the 
calculated solar models show a substantial differen- 
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Table 1 
Solar Models 



Input Parameters Results 





Mixing 




Overshoot 










"^crit 




CZ depth 


RUN 


Length 


Yo 


Zo/Xo 


(Hp) 


fmXk 


m/s) N= 


J, ^ 

JGSF 


M 


10"^s" 


M Lii/Uf ±0.0009 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(') 


/OA 

(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


BF 


1.841 


0.28167 


0.02768 


0.05 


1.0 


30 


1.5 


10 


1.00 


3.0 


199.2 


0.7096 


AM 


1.832 


0.28155 


02768 


0.05 


1.0 


30 


1.5 


10 


0.10 


3.0 


199.0 


0.7101 


CO 


1.827 


0.28155 


02768 


0.05 


1.0 


30 


1.5 


10 


0.01 


3.0 


200.9 


0.7109 


IF 


1.827 


0.28180 


02769 


0.02 


1.0 


30 


1.5 


10 


0.01 


3.0 


199.6 


0.7110 


HC 


1.827 


0.28175 


02769 


0.05 


1.0 


50 


1.5 


10 


0.01 


3.0 


200.6 


0.7100 


uu 


1.838 


0.28155 


0.02768 


0.05 


1.0 


30 


2.0 


10 


0.10 


3.0 


200.0 


0.7090 


DQ 


1.837 


0.28145 


0.02767 


0.05 


1.0 


50 


2.0 


10 


0.10 


3.0 


199.0 


0.7098 


RK 


1.828 


0.28185 


0.02769 


0.05 


1.0 


30 


2.0 


1.0 


0.10 


3.0 


202.1 


0.7097 


JF 


1.804 


0.28171 


02755 


02 


8 

U • (J 


30 


1.5 


10 


0.01 


3.0 


199.4 


0.7141 


QG 


1.800 


0.28199 


0.02756 


0.02 


0.8 


30 


1.5 


1.0 


0.01 


3.0 


199.4 


0.7137 


SF 


1.806 


0.28174 


0.02755 


0.02 


0.8 


30 


1.5 


10 


0.01 


1.5 


199.8 


0.7125 


TJ 


1.806 


0.28174 


0.02755 


0.02 


0.8 


30 


1.8 


10 


0.01 


1.5 


199.2 


0.7129 


VN 


1.806 


0.28207 


0.02757 


0.02 


0.8 


30 


1.5 


1.0 


0.01 


1.5 


200.2 


0.7126 


WE 


1.806 


0.28187 


0.02756 


0.02 


0.8 


30 


1.5 


1.0 


0.01 


1.5 


99.35 


0.7130 


EI 


1.725 


0.28125 


0.02686 


0.05 


0.0 


30 


1.5 


10 


0.10 


3.0 


200.1 


0.7241 


EC 


1.856 


0.28176 


0.02769 


0.05 


1.0 













4.747 


0.7089 


KB 


1.862 


0.28176 


0.02769 


0.02 


1.0 













3.648 


0.7078 


XB 


1.856 


0.28176 


0.02769 


0.25 


1.0 













202.4 


0.7094 


XC 


1.856 


0.28176 


0.02769 


0.24 


1.0 













148.5 


0.7094 


ND 


1.835 


0.28179 


0.02756 


0.02 


0.8 













3.358 


0.7117 


GD 


1.728 


0.28149 


0.02687 


0.05 


0.0 













2.999 


0.7246 


LA 


1.728 


0.28149 


0.02687 


0.02 


0.0 













2.461 


0.7246 


PA 


1.728 


0.28149 


0.02687 


0.00 


0.0 













2.151 


0.7246 



^Constant which multiplies the microscopic diffusion coefficients of Michaud & Proffitt 1993. 
'^Initial rotation velocity given to the models on the pre-main sequence. 
Power law index for the wind loss law. 
■^Constant which multiplies the estimate for the GSF diffusion coefficient. 

'^Constant which determines the efficiency of rotational mixing in the presence of mean molecular weight 
gradients. Small values of imply efficient mixing. 

^Saturation level in the angular momentum loss law. 

SDepth of the surface convection zone in units of Rq. In the Sun, the surface convection zone has a 
depth of 0.713 ± 0.003 Rq. 
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I't _ = 0.10, NO diffusion 




0.0 0.2 0.4 0.6 0.8 



Fig. 1. — Present day solar rotation curve for the 
rotating models EI (no diffusion), BF (/^ = 1.0), AM 
(/^ = 0.1) and CO (/^ = 0.01). 

tial rotation, with the core rotating many times faster 
than the surface. Below r = 0.2 Rq the rotation 
rate in model EI (no diffusion, = 0.1) rises steeply 
to match that of model AM (microscopic diffusion, 

= 0.1). Above r = 0.2 Rq, the rotation profile in 
model EI is nearly identical to model CO (/^ = 0.01). 
This clearly shows the inhibiting effects of microscopic 
diffusion on the rotational induced mixing in the outer 
layers of the Sun. Models with microscopic diffusion 
require a value of ten times smaller in order to 
achieve a similar rotation profile in the outer layers 
to the pure rotation models. 

Models with = 0.01 have the lowest core ro- 
tation rates, typically a factor of ~ 12 faster than 
the surface. Models with = 0.1 have core rota- 
tion rates ~ 25 times faster than the surface. Such 
high core rotation rates may be ruled out by the ob- 
served splittings of the lowest £ p-modes (Toutain & 
Frohlich 1992; Loudagh ei al. 1993) which suggest 
that the core of the Sun is rotating ~ 4 times faster 
than the surface. The lower core rotation rate in- 
ferred by Toutain & Frohlich and Loudagh may be a 
resolution problem, for even the lowest I modes are 
weighted such that they are ineffective at probing the 
rotation rate at a radius of ^ 0.2 Rq. At r = 0.2 i?0, 
models with = 0.01, 0.1 and 1.0 have rotation 
rates which are ~ 5, 8 and 23 times faster than the 
surface rate respectively. Thus, the = 1.0 case ap- 



pears to be ruled out on the basis of its extremely 
high core rotation rate. We note however, that if an 
angular momentum transport process which does not 
transport particles (such as magnetic fields, or waves) 
is active in the Sun, then it is not possible to rule out 
the = 1.0 case. The = 0.01 and = 0.1 are 
in reasonable agreement with the observations of the 
core rotation of the Sun. 

Observations of the £ = 10 - 60 yield information 
of the rotation velocity between r = 0.4 — 0.9 Rq)- 
Inversions of these observations indicate that the Sun 
rotates as a solid body over this range in radii (Lib- 
brecht k Morrow 1991, Goode ei al. 1991). All of 
our models contain a fairly substantial amount of dif- 
ferential rotation between the base of the convection 
zone and r = Q.A Rq. At r = 0.4 Rq, model CO 
(/^ = 0.01) has a rotation rate 2.9 times faster than 
the surface rate. Thus, none of our models agree 
with the rotation profile inferred by inversions of the 
£ = 10 - 60 p-modes. However, the inversion process 
is sensitive to observational errors, and the type of in- 
version procedure employed (see Goode et al. 1991). 
It is clear from observations of surface rotation rates 
in subgiants (Noyes et al. 1984; Pinsonneault et al. 
1989; Demarque & Guenther 1988) and in horizontal 
branch stars (Peterson 1985; Pinsonneault, Deliyan- 
nis & Demarque 1991) that a reservoir of angular 
momentum must be buried inside a star, for these 
stars rotate substantially faster than what can be ac- 
counted for by assuming a solid body rotation on the 
main sequence. For this reason we do not rule out any 
of our models based on their rotation profiles. How- 
ever, it appears that a more efficient angular momen- 
tum mixing mechanism is present in the Sun, which 
is not present in our models. Magnetic fields have 
been ignored in this study. Large magnetic fields in 
the radiative region of the Sun have the potential to 
efficiently mix angular momentum, and lead to a fiat 
solar rotation curve (see Charbonneau & MacGregor 
1993). 

The quadrupole moments and the associated solar 
oblateness have been calculated using the method of 
Fricke (1968, see Sweet 1950 for the original deriva- 
tion) modified to take into account the solar differ- 
ential rotation in latitude (Dicke 1970). Model BF 
(/^ = 1) has e = 2.8 x 10~®, which is 3(7 above the 
observation of Hill & Stebbins (1975), 2(7 above the 
observations of Dicke et al. (1985) and > 27(7 above 
the observations of Sofia et al. (1994). Thus, model 
BF is ruled out as a viable solar model on the basis of 
its high quadrupole moment. Model AM (/^ = 0.1) 
has e = 9.6 x 10~® which is in agreement with the 
measurements of Hill b, Stebbins (1975) and Sofia et 
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Fig. 2. — The radial distribution of the mass frac- 
tion (X) near the base of the convection zone in solar 
models which include microscopic diffusion (FC - no 
rotation; BF - = 1.0; AM - = 0.1; and CO - 
U = 0.01). 



Fig. 3. — The radial distribution of the Li divided by 
its initial value near the base of the convection zone in 
solar models which include microscopic diffusion (FC 
- no rotation; BF - = 1.0; AM - = 0.1; and 
CO - = 0.01). 



al. (1994). A similar statement holds for model CO 

iU = 0.01). 

It is clear from Figure 1 that models which include 
diffusion and a large value of show a distinct dis- 
continuity in uj at the base of the convection zone. 
This is due to the large gradient in /j, that develops 
due to the diffusion of '^He out of the convection zone. 
This /J, gradient effectively halts the rotational mixing 
in models with = 1.0 and 0.1. This may be seen 
in Figure 2 where we have plotted the mass fraction 
of (X) near the base of the convection zone for 
four models which include microscopic diffusion. In 
contrast to the high cases, model CO (/^ = 0.01) 
has a relatively small gradient in X at the base of 
the convection zone (r = 0.71 Rq)- This is due to 
the fact that the rotational mixing is able to smooth 
out some of the gradients composition caused by dif- 
fusion. Note that by r = 0.5 Rq, the hydrogen pro- 
file is very similar among the models. Figure 3 plots 
the radial profile of the ^Li abundance (divided by its 
initial abundance) for the same models. The rotat- 
ing models have been calibrated such that the surface 
^Li abundance is l/200th of its original value. The 
region of ^Li destruction is more extended in the ro- 
tating models than in the pure diffusion model due to 
the slow timescale of the rotational mixing. 



3.2. Effect of Rotation on the Microscopic 
Diffusion 

When rotational mixing is active, it tends to in- 
hibit the sorting effect of diffusion. This effect is best 
exemplified by examining the change in time of the 
surface '^He abundance in calibrated models (Figure 
4). It is best to concentrate on the surface abundance 
of ■^He because it is only infiuenced by particle trans- 
port; no nuclear destruction or creation contributes to 
its change. Figure 4 clearly demonstrates the interac- 
tion of rotational mixing with microscopic diffusion, 
and the effects of In model BF (with = 1.0) 
no rotational mixing occurs near the surface after 1.5 
Gyr and the rate of '^He depletion is virtually iden- 
tical to model FC (which has no rotational mixing). 
In contrast, in model AM (/^ = 0.1) the rate of '^He 
depletion does not match model FC until ~ 3.5 Gyr. 
The rate of '^He depletion is smallest (and relatively 
constant, at rate ~ 75% of the pure diffusion model) 
in model CO (/^ = 0.01). This suggests that in mod- 
els with small , the rotational mixing leads to an 
'effective' microscopic diffusion coefficient (used to de- 
rive Dm,i and 0^,2) which is about 30% smaller than 
the calculated coefficients near the surface of the star. 
Since the microscopic diffusion in the outer part of our 
model causes a change in the radius of the models, the 
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Fig. 4. — The lower panel shows the mass fraction of 
■*He (Y) at the surface of the Sun as a function of time 
for models which include diffusion (FC - no rotation; 
BF - = 1.0; AM - = 0.1; and CO - = 0.01). 
The upper panel shows the rate of change in Y . 

inhibiting effect of rotational mixing has important 
implications for effective temperatures of our stellar 
models. 

4. ^Li Depletion and Rotational History of 
the Solar Models 

The details of how the models arrive at their 
present state varies considerably. It is instructive to 
examine the variation in time of the surface rotation 
velocity and ^Li as these two quantities are readily 
observed in other stars. Figure 5 plots the surface 
rotation velocity (v) and ^Li depletion as a function 
of age for four models which include rotation. When 
examining this plot, it is important to remember that 
calibrated solar models have been calculated; the ro- 
tational mixing efficiency (/c) and wind constant { fk) 
have been adjusted to yield the observed surface rota- 
tion rate and ^Li abundance of the present day Sun. 

On the main sequence, the rate of decline in the 
surface v is steeper for models with large . This 
is because gradients in jj, inhibit the rotational mix- 
ing and decouple the core and the envelope. Thus, in 
models with large only the envelope spins down. 
This means that our calibrated models require a 
smaller wind constant, implying that the models with 
large spin up more during the pre-main sequence 




Age (Gyr) 

Fig. 5. — The surface rotation velocity (lower panel) 
and ^Li depletion as a function of time for models 
which include rotational mixing (EI - no diffusion; 
BF - = 1.0; AM - = 0.1; and CO - = 0.01) 

phase of evolution. Hence, beyond ~ 10 Myr, models 
with large have larger surface rotation velocities. 
It is interesting to note that the = 0.01 model 
has a nearly identical rotation velocity history to the 
= 0.10 pure rotation model. 

The effect of on the surface ^Li abundance 
is equally striking. In model BF (/^ = 1.0), the 
^Li abundance is virtually constant on the main se- 
quence, implying that very little rotational induced 
mixing occurs. However, observations of Li m the 
Pleiades and Hyades show that mam sequence Li de- 
pletion has occurred (Soderblom et al. 1993; Thor- 
burn et al. 1993). The = 0.10 and 0.01 models 
show significant main sequence Li depletion. Com- 
parisons to cluster observations are somewhat com- 
plex, as allowance must be made for a range in ini- 
tial rotation velocities. Most T-Tauri stars have ro- 
tation velocities near 10 km/s, though some have ro- 
tation velocities as high as 70 km/s (Bouvier et al. 
1993). To compare to the open cluster observations, 
Li destruction isochrones were calculated by evolving 
stellar models with masses of 0.6, 0.7, . . . 1.3 Mq and 
initial rotation velocities of 10 and 30 km/s. In all 
other respects, these models used parameters iden- 
tical to mo del BF (/^ = 1.0). The ^Li destruction 
isochrones are compared to ^Li observations in the 
Hyades (age 700 Myr) and M67 (age 4.5 Gyr) in Fig- 
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Fig. 6. — Li destruction isochrones for the = 1.0 
BF case are compared to the Hyades (lower panel) 
and M67 (upper panel) observations. Two sets of Li 
isochrones are shown - the solid line is for models 
with an initial rotation velocity of 10 km/s, while the 
dashed line represents the Li isochrone for a family of 
models with an initial rotation velocity of 30 km/s. 
It is clear that the 30 km/s isochrones are too low 
at the age of the Hyades (700 Myr), while the 10 
km/s isochrones are too high at the age of M67 (4 
Gyr). This is due to the fact that the models have 
little main sequence Li depletion. The observations 
are from Thorburn et al. (1993) for the Hyades, and 
from Hobbs & Pilachowski (1986) and Spite et al. 
(1987) for M67. 

ure 6. At the age of the Hyades, the 10 km/s Li 
destruction isochrones lie somewhat below the mean 
trend of the observations. The 30 km/s isochrones 
lie well below the mean trend. These models pre- 
dict a much larger dispersion (at a give Tgff ) than is 
observed in the Hyades. The lack of main-sequence 
^Li depletion causes the 10 km/s models to have much 
more ^Li than is observed in older clusters, such as 
M67. The lack of main-sequence ^Li depletion clearly 
rules out the = 1.0 case, implying that canonical 
estimates for the amount of inhibition in rotational 
mixing are too large. Reducing the inhibiting effect 
of mean molecular weights by a factor often, leads to 
noticeable main sequence ^Li depletion. 

Figure 7 illustrates the effect of changing the power 
law index N in the wind loss law on the models. As 
expected, models with higher N suffer significantly 




10'^ 10'' 10"' 10' 

Age (Gyr) 



Fig. 7. — Effect of changing the power law index N in 
the wind loss law on the surface V and ^Li abundance. 
Models AM {N = 1.5) and UU (N = 2) are shown. 
A Skumanich (1972) t~^/^ spin down law is shown in 
the bottom panel. 

more angular momentum loss. As a consequence, 
models with high N do not spin up as much, and 
so suffer less angular momentum loss on the main 
sequence. Thus, the surface ^Li depletion occurs 
somewhat earlier in models with high N. We have 
schematically illustrated a Skumanich (1972) t~^/^ 
spin down law on the Vrot plot. It is clear that model 
AM (N = 1.5, ffi = 0.1) spins down too rapidly on 
the main sequence. Thus, it is ruled out as a vi- 
able solar model. Model UU has the maximum N 
allowed (2) and follows the t~^/^ law reasonably well. 
In Figure 8 the sensitivity of the models to changing 
the saturation level in the angular momentum loss 
law is demonstrated. Models with low Wcrit spin up 
more doing the pre-main sequence and retain substan- 
tially faster rotation rates until an age of ~ 700 Myr. 
This is a natural consequence of lowering Wcrit, above 
which the power law dependence of the loss on uj is 
lowered. We note that past ~ 1 Gyr, the surface ro- 
tation velocities of these models are nearly identical. 
The ^Li depletion is somewhat different between the 
ages of ~ 60 — 200 Myr, which suggests that the ob- 
servations in the Pleiades (age ~ 70 Myr) may be able 
to distinguish between these models. 

In Figure 9 we show the consequence of varying the 
initial rotation velocity on the ^Li depletion and sur- 
face rotation velocity. For these calibrated models, a 
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Fig. 8. — Effect of clianging tlie saturation level in 
the angular momentum loss law on the surface V and 
^Li abundance. Models JF (wcnt = 3 x 10"^ s"!) SF 
= 1.5 X 10~® -s'^) ^re shown. 
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Fig. 9. — Effect of changing the initial rotation ve- 
locity on the surface velocity V , and ^Li abundance. 
Models CO (Vi = 30 km/s) and HC (Vi = 50 km/s) 
are shown. 
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Fig. 10. — Effect of changing the efficiency of the GSF 
instability in mixing material. Models JF (/gsf = 
10), and QG (/gsf = 1-0) are shown. 

higher initial velocity implies a larger wind constant 
as we require that the models rotate at 2 km/s at the 
solar age. During the early pre-main sequence phase 
(t ^ 25 Myr) the difference in surface rotation veloc- 
ity remains constant. However, once the star starts 
to spin down, the model with a higher initial rota- 
tion velocity loses angular momentum at a faster rate 
due to the higher wind constant. The ^Li depletion 
occurs somewhat earlier in models with higher ini- 
tial rotation velocities, though not by an measurable 
amount. 

Finally, Figure 10 shows the effect of modifying 
the GSF coefficients. Models with small /gsf do not 
spin up as much on the pre-main sequence and have 
lower surface rotation velocities through most of their 
pre-main sequence and early main sequence evolution 
(10 ^ 300 Myr). This is due to the fact that when 
the GSF instability is inefficient, the main sequence 
spin down of the star is slow, which requires that 
the star not spin up too much during the pre-main 
sequence (recall that we are examining solar models 
which have been calibrated to yield the surface solar 
rotation rate at the solar age). The depletion of ^Li is 
delayed when /gsf is lowered, due to the decreased 
efficiency of the GSF instability in mixing material. 
It is important to note that when the efficiency of the 
GSF instability is lowered, the mixing in higher mass 
stars is lowered relative to lower mass stars. 
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It is clear from Figures 7-10 that our mod- 
els display a complex behaviour in terms of surface 
^Li depletion and rotation velocities. However, a few 
general conclusions may be drawn. The power law 
index in the wind law (N) determines the rate of spin 
down on the main sequence. By varying it is pos- 
sible to achieve a Skumanich t~^/^ (1972) spin down 
on the main sequence. A low value for the angu- 
lar momentum saturation level leads to substantial 
higher pre-main sequence and ZAMS rotation veloc- 
ities. Models with different initial rotation velocities 
have nearly identical rotation velocities by an age of 
~ 150 Myr. 

5. Summary 

A number of solar models which include the com- 
bined effects of microscopic diffusion and rotationally 
induced mixing in the radiative regions of the model 
have been calculated. The models suffer angular mo- 
mentum loss from the surface. These are the first 
solar models which included all three of these effects. 
The models have been calibrated to match the ob- 
served solar radius, luminosity, surface Z/X, surface 
rotation velocity, and ^Li depletion. 

Helioseismology is able to determine the depth of 
the solar convection zone (Christensen-Dalsgaard ei 
al. 1991) to be 0.713 ±0.003 Rq. Solar models which 
include rotation and microscopic diffusion have con- 
vection zone depths of 0.710 Rq, providing a good 
match to the observed depth. Models which do not 
include diffusion have convection zones which are too 
shallow compared to observations. This has been 
found by other workers, who note that the convection 
zone depth in standard models may be made to agree 
with the observations by enhancing opacities near the 
base of the convection zone (Guzik & Cox 1993). 

Models which include microscopic diffusion build 
up gradients in the mean molecular weight near the 
base of the convection zone. Such gradients can in- 
hibit the rotationally induced mixing. The sensitivity 
of rotational mixing to gradients in the mean molec- 
ular weight is somewhat uncertain; we have parame- 
terized it by the non-dimensional multiplicative factor 

. A large implies a strong inhibiting factor due to 
gradients in mean molecular weight. The model with 

= 1.0 (the standard value, given by Kippenhahn 
1974) has very high internal rotation velocities which 
cause a large solar oblateness e = 2.8 x 10~®. Such a 
high value of the solar oblateness is ruled out by ob- 
servations (Sofia et al. 1994). Models with <^ 0.1 
are in agreement with the oblateness measurements. 
A further problem with assuming a high is that 



the model with = 1 depletes virtually no Li on 
the main sequence. This is in clear contrast with the 
observation of Li in the Pleiades and Hyades (Thor- 
burn et al. 1993; Soderblom et al. 1993). Thus, it is 
clear that gradients in the mean molecular weight in- 
hibit the rotational mixing by a smaller amount than 
previously thought. 

The rotational mixing inhibits the microscopic dif- 
fusion settling in the outer parts of the model. An 
examination of the surface '^He abundance (Figure 
4) reveals the amount of inhibition. Models which 
do not include rotational mixing require their micro- 
scopic diffusion coefficients (used to derive Dm,i and 
-Dm, 2) to be multiplied by 0.5 — 0.8 in order to obtain 
a similar depletion in the surface '^He as the models 
which include rotational mixing. Microscopic diffu- 
sion operates relatively unhindered in the core of the 
model. This may explain the findings of Chaboyer 
et al. (1992) that the amount of age reduction in 
globular clusters caused by diffusion depended on the 
method used. Isochrone fitting and A(B — V) meth- 
ods suggested age reduction of 2 - 3 Gyr, while the 
AV(TO — HE) method suggested an age reduction of 
0.5 - 1 Gyr. Both isochrone fitting and A(B — V) 
depend on the color (hence, radius) of the models, 
while AV(TO — HB) depends on the luminosity of the 
models. The amount of microscopic diffusion in the 
outer part of the star determines the change in radius 
(hence, color) in the models. The amount of diffusion 
in the core determines the effect of microscopic dif- 
fusion on the luminosity. Thus, the inhibiting effect 
of rotational mixing in the outer layers of the star is 
likely to reduce the age reduction found by isochrone 
fitting and A(B — V) by a factor of two while the age 
reduction found by AV(TO — HB) is unchanged. 

Observations of splittings in £ = 10 - 60 p-modes 
suggest that between the base of the convection zone 
and r = 0.4 Rq the Sun rotates as a solid body (Lib- 
brecht & Morrow 1991). Our model with = 0.01 
has the least amount of differential rotation of our 
models. At r = 0.4 i?0 it has an angular velocity 
2.9 times higher than the surface, which contradicts 
the solid body rotation inferred by inversions of the 
£ = 10 - 60 p-modes. This suggests that the trans- 
port of angular momentum is more efficient in the Sun 
than what we have modeled. We note however, that 
it is clear from the relatively rapid rotation of sub- 
giant stars (Pinsonneault et al. 1989) and of horizon- 
tal branch stars (Pinsonneault et al. 1991) that stars 
must have large reservoir of angular momentum which 
is transported to the surface once the stars leave the 
main sequence. All of the solar models have a core 
which rotates substantially faster than the surface. 
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The observations of the splitting in low I p-modes by 
Toutain & Frohlich (1992) and Loudagh et al. (1993) 
suggest that the rotation rate is ~ 4 times faster than 
the surface at r 0.2 Rq. At a similar depth, the 
= 0.1 and 0.01 cases have rotation rates 8 and 5 
times the surface rate respectively. Thus, they are in 
reasonable agreement with the splittings in the low I 
p-mode observations. 
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work. The comments of the anonymous referee im- 
proved the presentation of this paper. Research sup- 
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2136, NAGW-2469 and NAGW-2531. 
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